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Lagrangian perturbation theory for cosmological fluid describes structure formation in the quasi- 
nonlinear stage well. We present a third-order perturbative equation for Lagrangian perturbation 
with pressure in both the longitudinal and transverse modes. Then we derive the perturbative 
solution for simplest case. 
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CN ' I. INTRODUCTION 

X> 

^ ' The structure formation scenario based on gravitational instability has been studied for a long time. The Lagrangian 
perturbative method for the cosmological fluid describes the nonlinear evolution of density fluctuation rather well, 
t — ' Zel'dovich 0| proposed a linear Lagrangian approximation for dust fluid. This approximation is called the Zel'dovich 
approximation (ZA) J3, S El- El U- HI ZA describes the evolution of density fluctuation better than the Eulerian 
approximation LP, la B ■ After that, the second- and the third-order perturbative solution for dust fluid were derived 
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Recently the effect of the pressure in the cosmological fluid has been considered. At first, the effect of the pressure 
" is originated from velocity dispersion using the collisionless Boltzmann equation |l8lll9j| . Buchert and Dommguez [l8l | 
showed that when the velocity dispersion is regarded as small and isotropic it produces effective "pressure" or viscosity 
terms. Furthermore, they posited the relation between mass density p and pressure P, i.e., an "equation of state". 
Adlcr and Buchert have formulated the Lagrangian perturbation theory for a barotropic fluid. Morita and 
Tatekawa [2ll | and Tatekawa et al. |22t solved the Lagrangian perturbation equations for a polytropic fluid up to the 
second order. Hereafter, we call this model the "pressure model". 

Although the higher-order perturbative solution is expected to improve the approximation, there is a counterar- 
Q-f gument. Let us consider the evolution of the spherical void for dust fluid. Because the exact solution has already 
been derived, we can discuss the accuracy of the perturbative solutions. In this case, when we increase the order 
of Lagrangian approximation, contrary to expectation, the description becomes worse 0, H, Especially, when 
we stop the order of the perturbation until even order (the second-order), the perturbative solution describes the 
C$ • contraction of a void at a later. 

For the pressure model, although we do not know the exact solution for the evolution of the spherical void, the same 
problem may arise. In fact, according to a comparison between the first-order solution and the full-order numerical 
solution for the spherical symmetric model, the difference of these solutions obviously appears at a late time (density 
fluctuation S > 1 or 5 < -0.5) 

In this paper, we present the third-order perturbative equation for the pressure model. Then we derive the per- 
turbative solution for the simplest case, i.e., where the background is given by the Einstein-de Sitter (E-dS) Universe 
model and the polytropic index 7 = 4/3. Furthermore, we compare the evolution of the density fluctuation between 
the first-, second-, and third-order approximations for one-dimensional model. 

This paper is organized as follows. In Sec.^ we present a Lagrangian description of basic equations for cosmological 
fluids. In Sec. II I II we show perturbative equations and derive perturbative solutions for the pressure model up to 
the third order. In Sec. IIII Al we show the perturbative equations and solutions of the first- and the second-order 



perturbation. For the higher-order perturbation, we ignore the first-order transverse modes. In Sec. IIII Bl we present 
the perturbative equations for the third-order approximation. In general, it is extremely hard to solve the third-order 
perturbative equations. Therefore in Sec. IIII CI we derive the solutions for the simplest case (E-dS Universe model, 
7 = 4/3). 

In Sec. IIVI we introduce spacial symmetry. We consider a planar model and compare the evolution of the density 
fluctuation between the first-, second-, and third-order approximations. In Sec.0 we summarize our conclusions. 



'Electronic address: tatekawa@gravity.phys.waseda.ac.jp 



2 



In Appendix ^ we present the second-order perturbative equation for a spherical symmetric model in the pressure 
model. Because of mode-coupling of the first-order perturbations, it seems difficult to solve. 



II. BASIC EQUATIONS 

Here we briefly introduce the Lagrangian description for cosmological fluid. In the comoving coordinates, the basic 
equations for cosmological fluid are described as 

g + iv s >(i + *)} = o, (i) 

9v 1 , „ . a 1 _ 1 _ „ 

— + -(v-V x )v + -v = -g V X P, (2) 

at a a a ap 

V, x g = , (3) 

V x ■ g = -4irGp b a6 , (4) 

5^ P —^. (5) 
Pb 

In the Eulerian perturbation theory, the density fluctuation 5 is regarded as a perturbative quantity. On the other 
hand, in the Lagrangian perturbation theory, the displacement from homogeneous distribution is considered. 

x = q + s(q,t), (6) 

where x and q are the comoving Eulerian coordinates and the Lagrangian coordinates, respectively, s is the displace- 
ment vector that is regarded as a perturbative quantity. From Eq. 10, we can solve the continuous equation Eq. 
exactly. Then the density fluctuation is given in the formally exact form. 

J means the Jacobian of the coordinate transformation from Eulerian x to Lagrangian q. Therefore, when we derive 
the solutions for s, we can know the evolution of the density fluctuation. 
The peculiar velocity is given by 

v = as . (8) 

Then we introduce the Lagrangian time derivative: 

d I „ 

= +_„.V !B . 9 
at at a 

Taking divergence and rotation of Eq. @, we obtain the evolution equations for the Lagrangian displacement: 

V x ■ f s + 2-s~ K1P \ J-^V x j) = -AnGpbiJ- 1 - 1) , (10) 
\ a a z I 

V, x (s + 2^s) = 0, (If) 

where (') means the Lagrangian time derivative (Eq. ©)• To solve the perturbative equations, we decompose the 
Lagrangian perturbation to the longitudinal and transverse modes: 

s = V5 + s T , (12) 
V-s T = 0, (13) 

where V means the Lagrangian spacial derivative. 
Here we expand the Jacobian: 

J = 1 + Si t i + — (sijSjj — SijSj t i) + det (sjj) (14) 

= 1 + V 2 5 + \ {(V 2 S) 2 - S^Sji - sfjslt - 25,y4J 

+ det(S, ij + sl j ) . (15) 
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Because Eqs. (|10fl and (lllfl include the Eulerian spacial derivative, we change to the Lagrangian spacial derivative. 
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III. THE LAGRANGIAN PERTURB ATIVE SOLUTIONS 



A. The first- and second-order perturbative solutions 



From Eqs. I|1U|) and we obtain the first-order perturbative equations: 



V 2 [ §W + 2 ^SU - 4*G Pb SU - ^l£ll V 2 S (i)\ = Q (16) 
\ a ' a z 



V x (s T «+2-i T «) = 0. (17) 



The first-order solutions for the longitudinal mode depend on spacial scale. Therefore the solutions are described with 
a Lagrangian wavenumber. In this paper, we discuss only perturbative solutions in the E-dS Universe model. 

SM(K,t) = C + (K)D+{K,t) + C-{K)D-(K,t), (18) 



t- l / & J ±w {A\K\t-i+^Z) for 7 ^ 



D± ^ = ( t -i/«VSS7S=W far 7=f, (19) 



V^rG)^ = 4 1/3 

where J' means Bessel function. C ± (ii") is given by the initial condition. 
For the transverse mode, the solutions are same as for dust model: 

s T(1) oc const., r 1/3 . (20) 

The transverse mode does not have a growing solution. Therefore, in first-order approximation, the longitudinal mode 
dominate during evolution. Hereafter we consider only the longitudinal mode solutions for the first-order solutions. 

From Eqs. ifTUI) and ijTTJl . we obtain the second-order perturbative equations. For the longitudinal mode, the 
equation becomes 

V 2 f + 2*S^ - 4nGp b SW - IHELI^s^ ) = Q L(2) , (21) 



} L(2) 



2TrGp b 



o(l) c(l) 



V 2 S (r 



2 7 — 1 



k 7 ( 7 - l)p\ 



7-1 



v 2 v 2 s (1) v 2 s (1) ] - 



KIPI, 



7-1 



a \ /of 

For the transverse mode, after some arrangement, we can describe as follows 



c(l) c(l) 
°,ijk J ,ijk 



V 2 



,T(2) 



2 «.T(2) 



Q 



T(2) 



(22) 



(23) 
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Q T { 2) = ^7g_ (s^wt;. - sy/v-v-.sr - v^stv-s 1 ; y - <*\y.v-.s 



7 2 <?( 1 )v72 o(l) 



(l)n2(;(l) 

jk v ,ijfc 



(24) 



Here we notice for the second-order transverse mode solutions. In the pressure model, even if we consider only the 
longitudinal mode for the first-order, the second-order perturbation for the transverse mode appears. In dust model, it 
does not appear. Therefore, when we derive the third-order perturbative solutions, we must consider the second-order 
transverse mode. 

The second-order solutions are formally written as follows: 



5(2) 

? T(2) 



1 



\K\ 2 



1 



dt'G(K,t,t')Q L ^{K,t). 
dt'G T (t,t')Q T ^(K,t), 



\K\ 2 

G T {t 1 t l ) = 3(i' - i- 1 / 3 *' 4 / 3 ) . 
G L (K, t, t') depends on the equation of state. If 7 ^ 4/3 and v = 5/(8 — 67) is not an integer, we have 

G L (K, t, t') = --JL—(- 1+ ^j \~y^ 6 [j- v (A\K\t^^ 3 )j v (A\K\l/--r^ 3 ) 



(25) 

(26) 
(27) 



2 sin i/tt 



-J v (A\K\r^l^ J_ v (A\K\i!-~<+^ 



and if 7 = 4/3, 



G L (K, t, t') = --l^r'/V^r'Wt'W-iWc-w) , 



b(K) 



2b(K) 

— -B\K\*. 
36 1 1 



The second-order perturbative solution for the case of 7 = 4/3 in the E-dS universe model is described by 

S^(K) = JdK' Z L (K,K')E L (K,K',t), 
\K) = JdK'Zj(K,K')E T (K,K',t), 



(28) 

(29) 
(30) 

(31) 
(32) 



E L (K,K') = - 



{(K' ■ (K - K')) 2 - \K'\ 2 \K - K'\ 2 } 



Wi2 



(2tt) 3 \K\ 2 

+B (^\K'\ 2 \K - K'\ 2 (K' ■ (K - K 1 )) + \\K'\ A \K - K' 

+ ({K' ■ (K - K')f + 2\K - K'\ 2 (K' ■ (K - K 1 )) 2 ) }] . 



(33) 



B 



(2tt)3 \K\ 



K — K'\ (K 1 ■ (K - K')) 



[{K> ■ (K K'))K'j + \K- K'\ 2 K> \K'\ 2 (K, K]) (K' ■ (K K'))(K 3 K$] , (34) 



E L (K,K',t) = 



E T {K,K',t) 



where X)(e=±) means 



2b(K) 



E E 

;e=±) (®=±) 



C®(K')C®(K - K')t- 1 / 3(Bb< - K ')® b( - K - K ') 



(((Bb(K') ® b(K - K>) - \) 2 - b{K f) 



3E E 

(e=±) (®=±) 



C®{K')C®(K - K')t 



i\-t-l/3®b(K')®b(K-K') 



(®b(K') <g> b(K - K>) - l){®b(K') <g> b(K - K')) 



(e=±) 



(35) 
(36) 

(37) 
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B. The third-order perturbative equations 

The third-order perturbative equation becomes very complicated. Here we introduce scalar quantities generated by 
Lagrangian perturbation. 



Mi (A) = VA, 

Ma(A,B) = -((V-AXV-B)-^-^) 

jU 2 (A) ee /i2(A,A), 

/i 3 (A) = det(A (J ) , 



(38) 
(39) 

(40) 
(41) 



where A and _B are vector quantities. 

First, we show the longitudinal mode equation. As in the second-order perturbative equation, we separate the 
terms of the third-order perturbation from the others. Then the terms consisting of the first- and the second-order 
perturbations are collected to the source term Q L( ^\ 



.7-1 



(42) 



We consider the source term. Using Eqs. (|38 p -l]4i p . the terms are written as follows: 
Q L ^ = 4^4 Ml (5 ( i 1) ) 3 +/i 2 (^ 1) ,^ ) + s r (2) ) + M3(^ ) )" 



7-1 



(1) „(2) , T(2), 



(s;\s. 



•/4$(5 



-s 



+ 1 d j [nx(sf) d 3 (n x {sf) +l i2{Sf))} 
+d j {s%d k (fjL l {sf) + ^{sf))} 

+0, { (-71*1(5?) - 7Ma (Sj>) + ^f^M!(5«) 2 ) d m (S?) 
-S§dk {-7Mi(5j ) ) dj^P) - S§ 9^1(5?)} 

+d j [mis?) 5g a feM i(5j 3 )}" . 



(43) 



The transverse mode also seems complicated. However, if we neglect the first-order transverse mode, the evolution 
equation is described as 



(44) 



Q 



T(3) 



s\ 



(1) / S (2) 



a k ' 3 



S\ 



(1) / S (2) 



a J 



.fc 



-S ( 4 -w$ 



7-1 



*<S)} -{*( 



AirGpbS 



(i) 



K 7Pfc 1 / v 2o(l)o(l)c( 1 ) 



a- 



o(1)v72o(1) 



(45) 
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C. The third-order perturbative solutions — the simplest case 



As in the first- and the second-order solutions, the third-order solutions are described with the Lagrangian wavenum- 
ber. Following the method used in the second-order solutions, the third-order solution is given by this integration: 



= 

S-T(3) = 



\K\* 
1 



dt'G(K,t,t')Q L ^(K,t), 
dt'G T (t,t')Q T(3 \K,t), 



(46) 
(47) 



to 



Here we show the third-order perturbative solutions for simplest case, the case of 7 = 4/3 in the E-dS Universe 
model. In this case, the contribution of the gravitational terms and the pressure terms become identical: 



(48) 



Here we write these terms as M/a 3 (M — const.). From the terms described by the multiplication of and in 
Q L , the time evolution of a part of the third order perturbation is contributed a 

Ff(K,K",t) 



= J* dt'G(K,t,t') ^3 (sW(K",t')SW(K - K",t') 



BM 



Ab(K)b(K — K")\K\ 



E E E 

;e=±) (®=±) (0=±) 



r 1 ' 2 j dK'E L (K - K", K' - K") 

C®(K")C®(K' - K")C (K - K') 
2 63 - |) 2 - ^)((©6i ® 62 h - ±) 2 - bj) 



(49) 



61 = b{K"), b 2 = b(K' - K"), 63 = b(K - K'), 64 = b(K - K"), b 5 = b(K) . 



(50) 



On the other hand, from the terms described by the triplet of in Q L , time evolution of a part of the third order 
perturbation is contributed as: 

F 2 L (K,K\K",t) 

= J ' dt'G(K, t, t 1 ) ^3 (SW (K", t')S^ (K' - K" , t')S^ (K - K', t')) 

t -l/2 . 



M 

~2b(K)\K\ 2 

E E E 

e=±) (®=±) (o=±) 



C®(K»)C*(K' K")C®{K - gQ te6ieb20b3 
(861 e b 2 63 - i) 2 - b 2 



(51) 



When we consider the effect of pressure, even if s T « = 0, s T ( 2 ^ appears. Therefore, the contribution from the 
multiplication of and s T ^ also exists. The contribution from the multiplication of and s T ^ is given as 
follows: 

F£AK,K",t) 



f dt'G(K,t,t') ^3 (s^(K",t')^ 2 \K-K",t')) 

— j-i/2 /* d x' " T (K - K" K' - K") 

2b(K)\K\i j ■> 1 j 



E E E 

;©=±) (®=±) (0=±) 



C®(Jg ,/ )C (8 (fC / - K")C (K - K') 

(0& 2 h - i)(06 2 6 3 )((©6i & 2 &3 - ^) 2 - 6|) 



(52) 
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Next, we consider the third-order transverse mode s T ( 3 \ In the transverse mode, we must make the contribution 
not only from the multiplication of longitudinal modes, but also from the multiplication of and s T ^ 2 ^ : 

F?{K,K",t) 
<■* M 

M 



dt'G T (t,t>) ^3 (s^(K",t')S^(K-K'\t')^ 

r 1 ' 2 J dK 1 E L (K - K", K' - K") 

C®(K")C®{K - K')C (K' - K") 



2b(K- K")\K\ 



E E E 

(e=±) (®=±) (0=±) 



(®&i <8> 62 063 - \){®bi ®b 2 (Z>b 3 - ±)((0& 2 0&3 - \) 2 -b\)' 



F?(K,K',K",t) 

J' dt'G T (t,t') (^irV)^ 1 ^'-^^ 1 ^-^^ 



M 

l^ 1 



,-1/2 



E E E 

(e=±) (®=±) (0=±) 



(©61 b 2 63 - g)(©&i &2 h - i) 



fdt'CF{t,t') ^3 (&V{K»,1>)7% {2) {K-K"X)) 



1 ' - I dK' Ej(K-K",K' -K") 

C®{K")C®{K - K')C (K' - K") 



3M 

E E E 



(062 b 3 - i)(0& 2 &3)(0&i © b 2 63 - i)(06i © 62 63 - h) 



5b 1 ®b 2 0b 3 



(53) 



(54) 



(55) 



The third-order perturbative solutions for the case of 7 = 4/3 in the E-dS Universe model is described with Eqs. i|4P|) - 
(|54H . Even if we restrict ourselves to the simplest case, the third-order perturbative solutions become very complicated. 
For other case, although we can derive the solutions following similar procedures, the solutions may become difficult 
to analyze structure formation. 



IV. THE TIME EVOLUTION IN THE PLANAR MODELS 



In the previous section, we derived the third-order perturbative solutions for the simplest case. Although we 
restricted ourselves to the simplest case, the solutions are still complicated. For analyses of the perturbative solutions, 
we introduce spacial symmetry. If we consider a planar model, the nonlinear terms in the gravitational term disappear. 
The evolution equation for Lagrangian displacement becomes 



a a 2 (l + si)T +1 



(56) 



where s is Lagrangian displacement, not longitudinal mode potential pOj . From the expansion of Eq. (|56|l . we obtain 
perturbative equations. 

A „s 7-1 . , 

(57) 
(58) 



s 


+ 2^ (1) 
a 


- AirGpb^ - 


Kjpj 
a 2 


1 

,11 


= 0, 






■dP) 
s 


„ a -J 2 ) 
+ 2-s 
a 


- ATvGpt,^ - 


a 2 


,11 


kipT 

a 2 


1 

-(7- 




s 


+ 2^ (3) 
a 


- AnGpt s^ 3) - 


tvypV 

a 2 


^3) 
,11 


7-1 
KIPb 

a 2 


(7 4 





(59) 
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For the simplest case, i.e. the Einstein-de Sitter model and 7 = 4/3, using Eqs. (|49(l and i|51|l . we can describe the 
perturbative solution for a planar model: 



i* (1 W) 



C +{K)t- l / &+h ^ +C-(K)t- 1 / e - b W ) 
^ J dK' (K'(K-K') 2 )E L (K,K',t), 

^ j dK' f dK" [{K ■ K" ■ (K - K")){K" - K')K'\ 



(60) 
(61) 



E E E 

(®=±) (®=±) (0=±) 

140 
"~2T 



C®{K")C®(K' - K")C (K - K') 
((®&2 h - ±) 2 - &!)((©&! ®6 2 06 3 - ±) 2 - 6|) 1 

J dK' J dK" [K'K"(K - K' - K") 2 ] 

1 (K- K'} 



3bl®&20&3 



E E E 

(®=±) c®=±) (0=± 



C®(X")C«(K' - X")C«(K - K')^ biQb20b3 



(®&i e 6 2 6 3 - |) 2 - &§ 



(62) 



Here we compute power spectra V(k) = (5 2 (k)J of density fluctuation. We choose the initial spectrum index as 

n = +1,0, and —1. The initial amplitude of density fluctuation is set as <5max — 10 -3 at a = 1. The Jeans 
wavenumber is given by hand, Kj = 80. In general, the Jeans wavenumber depends on time. However, the Jeans 
wavenumber will still be constant in our calculations, because we choose the polytropic index 7 = 4/3. 

Initial conditions for two independent quantities are required to determine . To determine , we set the initial 
density fluctuation 5o and the initial peculiar velocity "Uo by growing mode solution in ZA. The procedure was shown 
by Morita and Tatekawa [2l| . 

In ZA, the density fluctuation diverges at a ~ 1000. Here we observe the difference between first-, second-, and 
third-order approximations. Therefore, we calculate the evolution just before shell-crossing, i.e., density divergence. 

Fig-G]shows the power spectra for n — +1, 0, —1 at a = 8000 (n = 1), a = 3000 (n = 0), and a = 1000 (n = —1). In 
these figures, we take an ensemble average of 900 samples. For n = 1 and n = 0, a difference between the results by 
the first-, second-, and third-order Lagrangian approximations appears. The effect of the higher-order approximation 
suppresses the evolution of the density fluctuation. However, for n = — 1, because the initial power in the large 
wavenumber component is small, the pressure does not surpress the evolution well. Therefore, the difference between 
the Lagrangian approximations still be very small even just before shell-crossing. 

In the one-dimensional model, the pressure only affects the nonlinear effect. However, in the three-dimensional 
realistic model, the gravitational force also affects the nonlinear effect, and the difference between first-, second- 
, and third-order approximation obviously appears. In fact, according to the comparison between the first-order 
approximation and full-order numerical calculation, the difference becomes large in the strongly nonlinear region [24| . 



V. SUMMARY AND CONCLUDING REMARKS 



In this paper, we showed the third-order Lagrangian perturbative equations for the cosmological fluid with pressure. 
Then we derived third-order perturbative solutions for a simple case. 

In our analysis for a planar model, the effect of the third-order perturbation seems very small even at the nonlinear 
stage. However, our result does not show that we can ignore the third-order perturbation easily. As mentioned 
in Sec. IIVI when we introduce planar symmetry, the nonlinear term of the gravitational force disappears. When we 
consider the effect of nonlinear pressure and gravitational force, the third-order perturbation is expected as a powerful 
tool to treat high-density regions. 

Recently several dark matter models have been proposed [2q| . If the interaction in some kind of dark matter can be 
described by the effective pressure, we can examine the behavior of the density fluctuation in a quasi-nonlinear stage. 
Furthermore, when we compare the observations and the structure that is formed by using the pressure model, we can 
delimit the nature of the dark matter. Especially when we consider high-density regions, the third-order perturbative 
solutions may become useful. 



9 



Acknowledgments 

We are grateful to Kei-ichi Maeda for his continuous encouragement. We would like to thank Hajime Sotani for 
useful discussion and comments regarding this work. The numerical calculation was in part carried out on the general 
common-use computer system at the Astronomical Data Analysis Center, ADAC, of the National Astronomical 
Observatory of Japan and Yukawa Institute Computer Facility. This work was supported by a Grant-in-Aid for 
Scientific Research Fund of the Ministry of Education, Culture, Sports, Science, and Technology (Young Scientists 
(B) 16740152). 



[1] Ya. B. Zel'dovich, Astron. Astrophys. 5, 84 (1970). 

[2] V. I. Arnol'd, S. F. Shandarin, and Ya. B. Zel'dovich, Geophys. Astrophys. Fluid Dynamics 20, 111 (1982). 
[3] S. F. Shandarin and Ya. B. Zel'dovich, Rev. Mod. Phys. 61, 185 (1989). 
[4] T. Buchert, Astron. Astrophys. 223, 9 (1989). 

[5] P. Coles and F. Lucchin, Cosmology: The Origin and Evolution of Cosmic Structure (John Wiley & Sons, Chichester, 
1995). 

[6] V. Sahni and P. Coles, Phys. Rep. 262, 1 (1995). 

[7] D. Munshi, V. Sahni, and A. A. Starobinsky, Astrophys. J. 436, 517 (1994). 
[8] V. Sahni and S. F. Shandarin, Mon. Not. R. Astron. Soc. 282, 641 (1996). 
[9] A. Yoshisato, T. Matsubara, and M. Morikawa, Astrophys. J. 498, 48 (1998). 
[10] J. D. Barrow and P. Saich, Class. Quantum Grav. 10, 79 (1993). 

[11] F. R. Bouchet, R. Juszkiewicz, S. Colombi, and R. Pellat, Astrophys. J. 394, L5 (1992). 

[12] F. R. Bouchet, S. Colombi, E. Hivon, and R. Juszkiewicz, Astron. Astrophys. 296, 575 (1995). 

[13] T. Buchert, Mon. Not. R. Astron. Soc. 254, 729 (1992). 

[14] T. Buchert and J. Ehlers, Mon. Not. R. Astron. Soc. 264, 375 (1993). 

[15] T. Buchert, Mon. Not. R. Astron. Soc. 267, 811 (1994). 

[16] P. Catelan, Mon. Not. R. Astron. Soc. 276, 115 (1995). 

[17] M. Sasaki and M. Kasai, Prog. Theor. Phys. 99, 585 (1998). 

[18] T. Buchert and A. Dommguez, Astron. Astrophys. 335, 395 (1998). 

[19] J. Binney and S. Tremaine, Galactic Dynamics (Princeton University Press, Princeton, NJ, 1987). 

[20] S. Adler and T. Buchert, Astron. Astrophys. 343, 317 (1999). 

[21] M. Morita and T. Tatekawa, Mon. Not. R. Astron. Soc. 328, 815 (2001). 

[22] T. Tatekawa, M. Suda, K. Maeda, M. Morita, and H. Anzai, Phys. Rev. D66, 064014 (2002). 

[23] T. Tatekawa, astro-ph/0412025 

[24] T. Tatekawa, Phys. Rev. D70, 064010 (2004). 

[25] J. P. Ostriker and P. Steinhardt, Science 300, 1909 (2003). 



APPENDIX A: SECOND-ORDER EQUATION FOR SPHERICAL SYMMETRIC MODEL WITH 

PRESSURE 



When we consider spherical collapse or expansion of a spherical void, we introduce spherical symmetry. 
The second order perturbative equation (Eq. 1)21(1 ) is changed as follows: 



V-' ! C + 2^C - 4ttGM - ^(p h )V 2 C 
a a z dp 



2nGp b 
dP 



-2S'{S' + 2rS") 



dp 



iPb) 



(-25' + 2rS" + r 2 S"') 2 



6(S') 2 - 12rS'S" + 6r 2 (S") 2 + r 4 (S'") 2 



S'(-4S' + 8rS" + 2r 2 S 



2C"\ i „2c« 



S"(-4S" + 2rS"' + r 2 S""^ 



d 2 P 
dp 2 



(Pb)pb 



{AS" 1 + rS"")(2S' + rS") (-25' + 2rS" + r 2 S'") 2 



10 



. „ S'(S / + 2rS") 
-AnGpb— — — 2 

1 H P 

—=r=-(Pb) [2(5') 2 - ±rS'S" + 2r 2 {S") 2 + 8r 3 S"S"' + 2r\{S"'f + S" S""j\ 
r 4 dp 

J >VP (p b ) Pb [4(5') 2 - SrS'S" + 4r 2 ((S") 2 + S'S'") + 2r 3 (S'S"" + 4S"S"') + r\S"S"" + (S'") 2 )] , (Al) 



,■4 



dp 2 



where S and £ are the first- and the second-order Lagrangian perturbation, respectively. Because the mode-coupling 
of the first-order perturbation becomes very complicated, it seems very difficult to solve the equation IjAljl . 
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FIG. 1: Power spectra of density fluctuation computed from Lagrangian first-, second-, and third-order approximations, (a) 
The initial spectrum index is n = 1. The power spectra are at a — 8000. The difference between the results by the first-, 
second-, and third-order Lagrangian approximation appears, (b) The power spectra are at a — 3000 for n = 0. In this case, 
the effect of the third-order approximation also appears, (c) The power spectra are at a — 1000 for n = — 1. The difference 
from Lagrangian approximation will still be very small just before shell-crossing. 



